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The effect of interchain interaction on the pairing symmetry competition in organic 
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We investigate the effect of interchain repulsive interaction on the pairing symmetry competi- 
tion in quasi-one-dimensional organic superconductors (TMTSF)2X by applying random phase 
approximation and quantum Monte Carlo calculation to an extended Hubbard model. We find 
that interchain repulsive interaction enhances the 2&f charge fluctuations, thereby making the 
possibility of spin-triplet /-wave pairing dominating over spin-singlet d-wave pairing realistic. 

KEYWORDS: (TMTSF)2X, superconductivity, spin-triplet pairing, interchain interaction, charge and spin 
fluctuations 



Possible occurrence of unconventional superconduc- 
tivity in organic conductors has been of great inter- 
est recently. Microscopically understanding the mecha- 
nism of pairing in those materials is an intriguing theo- 
retical challenge. Among the various candidates of un- 
conventional superconductors, the pairing mechanism 
of quasi-one-dimensional (qlD) organic superconductors 
(TMTSF) 2 X {X = PF 6 , C10 4 , etc.), so called the 
Bechgaard salts, 1 ' 2 has been quite puzzling. Namely, 
since superconductivity lies right next to the 2fcp spin 
density wave (SDW) phase in the pressure-temperature 
phase diagram, a spin-singlet ri-wave-like pairing (shown 
schematically in Fig. 1(a)) is expected to take place as 
suggested by several authors. 3-5 However, an unchanged 
Knight shift across T c 6 and a large H C 2 exceeding the 
Pauli limit suggest a realization of spin-triplet pairing. 
As for the orbital part of the order parameter, there have 
been NMR experiments suggesting the existence of nodes 
and thus unconventional pairing, 8 although a thermal 
conductivity measurement suggests absence of nodes for 
(TMTSF) 2 C10 4 . 9 

As a possible solution for this puzzle of spin-triplet 
pairing, one of the present authors has phenomelogically 
proposed that triplet /-wave-like pairing (whose gap is 
shown schematically in Fig. 1(b)) may take place due to 
a combination of quasi-lD (disconnected) Fermi surface 
and the coexistence of 2fcp spin and 2fcp charge fluctua- 
tions. 10 Namely, due to the disconnectivity of the Fermi 
surface, the number of gap nodes that intersect the Fermi 
surface is the same between d and /. Moreover, if the 2fc^ 
spin and charge fluctuations have about the same mag- 
nitude, spin-singlet and spin-triplet pairing interactions 
have close absolute values (with opposite signs) as will 
be explained later. In such a case, spin-triplet /-wave 
pairing should be closely competitive against singlet ci- 
wave pairing. As for other possibilities of triplet pairing, 
the p-wave state in which the nodes of the gap (Fig. 1(c)) 
do not intersect the Fermi surface has been considered 
from the early days, 11-13 but from a microscopic point of 
view, spin-triplet pairing interaction has a negative sign 
for the momentum transfer of 2kp unless spin fluctua- 
tions are highly anisotropic, so that a gap that changes 
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Fig. 1. Candidates for the gap function of (TMTSF) 2 X are 
schematically shown along with the Fermi surface (solid curves). 
The dashed lines represent the nodes of the gap, whose ki, depen- 
dence is omitted for simplicity. (For the actual ki, dependence, 
see Fig. 3.) We call the gap in fig.(a)((b)) d-wave (/-wave) in the 

sense that the gap changes sign as H 1 — (H 1 1 — ) along 

the Fermi surface. 



sign between the left and right portions of the Fermi sur- 
face is unlikely to take place. 14 A similar phenomelogi- 
cal proposal of /-wave pairing in (TMTSF^X has also 
been given by Fuseya et al. 15 Experimentally, the /-wave 
scenario due to the coexistence of 2/cf spin and charge 
fluctuations is indirectly supported by the observation 
that 2kp charge density wave (CDW) actually coexists 
with 2Uf SDW in the insulating phase lying next to the 
superconducting phase. 16 ' 17 

As for microscopic theories for the pairing competition, 
we have previously shown using a ground state quan- 
tum Monte Carlo method that /-wave strongly domi- 
nates over p-wave in the Hubbard model that considers 
only the on-site repulsive interaction. 18 More recently, 
we have shown, by applying random phase approxima- 
tion (RPA) to an extended Hubbard model, that /-wave 
pairing can indeed dominate over ri-wave pairing when 
we have large enough second nearest neighbor repulsion 
V', 19 which has been known for some years to have the 
effect of stabilizing 2Uf CDW configuration. 20,21 To be 
more precise, the condition for /-wave dominating over 
d-wave is to have V ~ U/2 (where U is the on-site repul- 
sion) or larger V because 2kp spin and 2Uf charge fluc- 
tuations have the same magnitude for V' — U/2 within 
RPA. A similar condition for /-wave being competitive 
against d-wave has also been obtained in a recent renor- 
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Fig. 2. The model of the present study is shown. 



malization group study. 22 Although these results do sug- 
gest that /-wave pairing can indeed be realized in mi- 
croscopic models, the condition that the second nearest 
neighbor repulsion being nearly equal to or larger than 
half the on-site repulsion may not be realized so easily 
in actual materials. In the present study, we consider 
a model where the interchain repulsion is taken into ac- 
count, which turns out to give a more realizable condition 
for /-wave dominating over cZ-wave due to the enhance- 
ment of 2Uf charge fluctuations. After completing the 
major part of this study, we came to notice that a sim- 
ilar conclusion has been reached quite recently using a 
renormalization group approach. 23 

The model considered in the present study is shown in 
Fig. 2. In standard notations, the Hamiltonian is given as 



H 



E 

<i,j>,a 



tij c ia c j<? 



+ u z2 n n n n + Y " 



<i,j> 



where c\ a creates a hole (note that (TMTSF^X is ac- 
tually a 3/4 filling system in the electron picture) with 
spin a =], I at site i. As for the kinetic energy terms, we 
consider nearest neighbor hoppings = t in the (most 
conductive) a-direction and Uj = t± in the ^-direction. 
t is taken as the unit of energy, and we adopt t± = 0.2t 
throughout the study. U and Vij are the on-site and the 
off-site repulsive interactions, respectively, where we take 
into account the nearest neighbor interchain repulsion V± 
in addition to the intrachain on-site (U), nearest (V), 
next nearest (V), and third nearest (V") neighbor re- 
pulsions considered in our previous study. 19 

Within RPA, 24-26 the effective pairing interactions for 
the singlet and triplet channels due to spin and charge 
fluctuations are given as 



V s (q) = U + V(q) + -U 2 Xs (q) 
- l -(U + 2V(q))\ c {q) 



V t {q) = V{q)--U\ s {q) 
-\{U + 2V{q)f Xc {q), 



(1) 



(2) 



where 



V(q) = 2V<^q x +2V'cos(2q x )+2V''cos(3q x )+2V±cos{q y 

(3) 



Here, Xs and Xc are the spin and charge susceptibilities, 
respectively, which are given as 

Xo(q) 



Xs(q) 



Xc{q) 



xo(q) 



l + (U + 2V(q))xo(q)' 
Here xo is the bare susceptibility given by 

1 ^ /(ep + q )~ f(e p ) 

xo(q) = ^l^ — — , - 



(4) 



p e P - e p + q 



with e^, = — 2tcosk a — 2t j_ cos fc& — /i and f(cp) — 
l/(exp(ep/T) + 1). xo peaks at the nesting vector Q 2kp 
(= (7r/2,7r) here) of the Fermi surface. 

To obtain T c , we solve the linearized gap equation 
within the weak-coupling theory, 



A s '*A s '*(fc) 



Ty*(fc-y) tmh ^ /2) A^ 



(5) 

The eigenfunction A s,t of this eigenvalue equation is the 
gap function. The transition temperature T c is deter- 
mined as the temperature where the eigenvalue A reaches 
unity. Note that the main contribution to the summation 
in the right hand side comes from k — k' ~ Q 2kp be- 
cause V s,t (q) peaks around q = Q 2kp . Although RPA is 
quantitatively insufficient for discussing the absolute val- 
ues of T c , we expect this approach to be valid for studying 
the competition between different pairing symmetries. 

Now, from eqs.(3) and (4), it can be seen that 
Xc{Q 2kF ) = Xs(Q 2kF ) holds when V' + V± = U/2, which 
in the absence V± of course reduces to the condition 
V' = U/2 obtained in our previous study. This in turn 
results in V s {Q 2kp ) — —V t (Q 2kF ) for the pairing inter- 
actions apart from the first order terms as can be seen 
from eqs.(l) and (2). Thus, considering the fact that the 
number of nodes intersecting the Fermi surface is the 
same between d and /, the condition for /-wave being 
competitive against d-wave should be V + V± ~ U/2. 
The possibility of this condition being satisfied in actual 
materials is realistic since V± can be comparable with the 
intrachain off-site repulsions due to the fact that the lat- 
tice constant in the a— and b— directions are of the same 
order. An intuitive picture here is that V± tends to "lock" 
more firmly the 2kp charge configuration induced by V', 
so that 2kp charge fluctuations are enhanced, thereby 
stabilizing the spin-triplet /-wave state. 

Bearing the above analysis in mind, we now move on 
to the RPA calculation results for the pairing symmetry 
competition between /- and d-waves. We first focus on 
the case where the parameter values satisfy the condition 
for Xc(Q 2kF ) = Xs{Q 2kF ), that is when V' + V± = U/2 
holds. Here we take U = 1.7, V = 0.8, V = 0.45, 
V" = 0.2, and V± = 0.4 in units of t. Note that V 
is much smaller than U/2. As expected, the singlet pair- 
ing having the largest eigenvalue A has a d-wave gap, 
while the triplet pairing with the largest A has a /- 
wave gap, as seen in Fig. 3. In Fig. 4, we plot A as func- 
tions of temperature for cZ-wave and /-wave pairings. 
The two pairings closely compete with each other, but 
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(a) triplet 



(b) singlet 




Fig. 3. The gap functions having the largest eigenvalue for the 
(a)triplet and (b)singlet pairing channels. The parameter values 
are taken as U = 1.7, V = 0.8, V = 0.45, V" = 0.2, V ± = 0.4 
and T = 0.011 (=T C of the /-wave pairing). The dark dashed 
curves represent the nodes of the gap, while a pair of light dotted 
curves near k a = ±tt/4 is the Fermi surface. 
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Fig. 5. T c plotted as a function of V± for U = 1.7, V = 0.8, 
V' = 0.45, and V" = 0.2. The solid (dashed) curve represent the 
/ (d)-wave regime. 




0.4 I , , , , T^J 

0.01 0.012 0.014 0.016 0.018 0.02 

T 

Fig. 4. The largest eigenvalue in the singlet and the triplet chan- 
nels are plotted as functions of temperature for U = 1.7, V = 0.8, 
V = 0.45, V" = 0.2, and V± = 0.4. 

/-wave pairing dominates over d-wave pairing and gives 
a higher T c . /-wave not being degerate with d-wave even 
for V' + V± = U/2 is due to the effect of the first or- 
der terms in eqs.(l) and (2) as discussed in our previous 
study. 19 

To look into the effect of V± on the /- vs. d- com- 
petition in more detail, we plot T c along with the pair- 
ing symmetry as a function of V±_ in Fig. 5. The pair- 
ing symmetry is /-wave and T c increases with V±_ for 
V± > U/2 — V'(—0A here), while the pairing occurs 
in the <i-wave channel with a nearly constant T c for 
V± < U/2 — V . The increase of the /-wave T c is due 
to the enhancement of 2kp charge fluctuations with in- 
creasing Vj_. 

Finally, in order to check the validity of RPA, we 
have performed auxiliary field quantum Monte Carlo 
(AFQMC) calculation 27-29 for the same extended Hub- 
bard model. Let us first briefly summarize this method. 
In AFQMC, the density operator is decomposed into 
the kinetic energy part and the interaction part using 
Trotter-Suzuki decomposition, 30,31 and we perform dis- 
crete Hubbard-Stratonovich transformation 27 for the in- 



teraction part. The summation over the Stratonovich 
variables are taken by Monte-Carlo importance sam- 
pling. Using this method, correlation functions and sus- 
ceptibilities can be calculated for finite size systems (16 
sites in the a-direction and 4 sites in the 6-direction=64 
sites in the present study), and the results are exact 
within the statistical errors. A defect of this approach 
is that we cannot go down to very low temperatures 
in the presence of off-site repulsions such as V, V and 
V± due to the negative sign problem, so that it is diffi- 
cult to look into the pairing symmetry competition itself. 
Nevertheless, we can check the validity of RPA at mod- 
erate temperatures of the order of Q.lt. Here we com- 
pare the values of Xs(Q2k F ) an d Xc{Q2k F ) calculated 
by AFQMC at T = 0.25,fixing V = 0.9, V" = and 
Vi = 0.3. In Fig. 6, we show the "phase diagram" in 
U — V plane, where we find that the AFQMC bound- 
ary for Xc{Q 2 k F ) = Xs(Q 2 k F ) is very close to the RPA 
boundary V' + V± = U/2. This result suggests that the 
RPA condition for Xc{Q2k F ) = Xs{Q2k F ) IS reliable at 
least at moderate temperatures. 

To summarize, we have studied the pairing symmetry 
competition in a model for (TMTSF^X which considers 
not only the intrachain repulsions but also the interchain 
repulsion. We find that the possibility of satisfying the 
condition for realizing /-wave pairing becomes more real- 
istic in the presence of the interchain repulsion. It would 
be an interesting future study to investigate whether this 
condition is actually satisfied in (TMTSF^X using first 
principles or quantum chemical calculations. Experimen- 
tally, it would be interesting to further confirm spin- 
triplet pairing by using probes complementary to those 
in the previous studies, 6,7 for example, a phase sensitive 
tunneling spectroscopy study 32 with 33 or without 3 ' 35 
applying a magnetic field, or those based on a newly 
developed theory for triplet superconductors, which has 
been proposed by one of the present authors. 36 
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Kobayashi for motivating us to study the effect of in- 
terchain repulsion. He also thanks J. Suzumura and Y. 
Fuseya for valuable discussion. Part of the numerical cal- 
culation has been performed at the facilities of the Su- 
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Fig. 6. AFQMC result for the competition between Xs(Q2k F ) 
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